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THE  ACCURACY  OF  CURVE  FITTING  PROCEDURES 


FOR  THE  DETERMINATION  OF  THE  RATE  CONSTANT 
FOR  REACTIVATION  OF  AN  ACYLATED  ESTERASE 


INTRODUCTION  | 

; 

j 

The  solution  of  equations  encountered  In  enzyme  kinetics  often 
requires  the  use  of  curve  fitting  procedures  CID.  Recent  work  has 
emphasised  that  rigorous  testing  of  such  procedures  Is  necessary  before 
they  can  be  utilised  confidently  In  an  experimental  situation  C2-4D.  One 
example  of  an  application  In  which  curve  fitting  procedures  fall  to  give 
reliable  results  Is  given  In  the  present  paper.  The  case  under  considera- 
tion Is  that  of  an  esterase  which  has  been  Inactivated  by  an  acylatlng 
Inhibitor.  On  removal  of  free  Inhibitor  the  enzyme  recovers  activity 
according  to  equation  (1). 

y = A-Be"'^*^  (1) 

where  y Is  the  enzymic  activity  at  time  t and  A,  B and  k are  constants  C5D. 

Specifically,  A Is  the  activity  of  completely  re-actlvated  (uninhibited) 
enzyme  and  (A-B)  Is  the  activity  at  t ~ 0.  If  A Is  known,  equation  (1) 

can  be  reduced  to  a linear  form  (equation  2)  and  the  rate  constant  k and 
the  constant  B determined  by  least-squares  linear  regression  on 
In(A-y)  vs  t. 


In(A-y)  = InB-kt  (2) 

Often  however  the  value  of  A Is  not  known  accurately  and  k can  then  be 
determined  only  by  more  complex  curve  fitting  procedures  1161].  A summary  of 
the  procedures  likely  to  be  adopted,  and  their  accuracy.  Is  given  below. 
Simulated  and  experimental  data  were  examined.  The  experiments  showed  that 
If  all  three  parameters  were  determined  by  various  curve  fitting  methods, 
very  large  errors  could  occur.  Even  In  cases  where  apparently  good  results 
were  obtained,  error  analyses  showed  that  confidence  ranges  for  k were  very 
wide.  If  A was  taken  as  a known  constant,  however,  the  fitted  values  of 
the  other  two  parameters  were  much  closer  to  their  correct  values  and  the 
confidence  ranges  for  k were  much  narrower. 


METHODS 


i 

"! 

i 


Derivation  of  Simulated  Data 


A set  of  random  nuinbers  between  0 and  1 was  generated  by  a PDF  10 
computer.  Algebraic  manipulation  of  the  numbers  produced  two  random 
error  sets;  one  In  which  the  numbers  varied  about  zero  (set  El)  and  the 
other  In  which  the  numbers  varied  about  1 (set  E2).  A third  random  number 
set  (E3)  was  derived  from  a different  set  of  computer-generated  random 
numbers  between  0 and  1.  The  three  error  sets  are  listed  In  Table  1. 
Appropriate  values  were  assigned  to  the  constants  A,  B and  k In  equation 
(1),  and  8 theoretical  values  of  y thereby  obtained  (Table  II).  The 
number  of  theoretical  values  (8)  was  chosen  as  being  a likely  number  of 
observations  In  an  experimental  situation.  The  simulated  experimental 
values  of  y were  then  obtained  from  the  theoretical  values  using  one  of 
the  error  sets  In  Table  I.  These  simulated  values  and  their  method  of 
calculation  are  also  given  In  Table  II.  Several  sets  of  simulated  data 
were  obtained  In  this  manner  (Table  II). 


Experimental  Data 


Bovine  erythrocyte  acetylcholinesterase  (acetylcholine  hydrolase, 

EC  3. 1.1.7)  was  Inhibited  with  neostigmine  bromide  (a  phenolic  ester  of 
dlmethylcarbamlc  acid) . The  solution  of  dime thy Icarbamylated  enzyme  was 
diluted  extensively  to  reduce  the  concentration  of  neostigmine  to  a non- 
Inhlbltory  level  and  the  enzymic  activity  was  measured  at  various  times 
after  dilution  by  the  method  of  Ellman  et  al.  173.  Full  experimental 
details  will  be  published  later  (R.M.  Dawson,  manuscript  In  preparation). 

In  the  presence  of  1.8  mM  choline  chloride,  the  half-life  for  recovery  of 
enzymic  activity  was  observed  to  be  approx.  13  min,  and  the  activity  at 
180  min  (14  half-lives)  was  taken  to  represent  100%  enzymic  activity 
(constant  A In  equation  1).  The  experimental  data  (Table  III)  were  plotted 
In  accordance  with  equation  (2)  and  the  constants  B and  k thereby  deter- 
mined by  least-squares  linear  regression.  Calculated  values  of  y were 
then  determined  for  each  value  of  t under  study  using  equation  (1)  and  the 
"known"  values  of  A,  B and  k,  and  the  calculated  values  compared  with  the 
experimental  values  (Table  III). 


Procedures  for  Curve  Fitting 

All  calculations  were  done  cn  the  computer  using  programs  written  In 
BASIC  or  FORTRAN.  Six  methods  of  curve  fitting  were  Investigated  In  order 
to  find  best-fit  values  of  A,  B and  k.  They  were  : 


Minimise  the  standard  error  of  the  least-squares 
straight  line  through  the  data  points  (ln(A-y),t)  C83. 

Minimise  an  approximation  to  the  standard  error. 

« 2 

Minimise  the  sum  of  squares  I(y-y)  C13. 

Minimise  ZC(y-y)^/y3. 


Method  1 : 

Method  2 : 
Method  3 : 
Method  4 : 


1 


i 


2 


FT 


Method  5 : Minimise  EC (y-y)^/(y+y)^D. 

Method  6 ; Minimise  E | y-y | . 

In  the  above,  y are  the  "experimental"  values  at  known  times  t 
(Tables  II  and  III)  and  y the  values  calculated  from  equation  (1).  Two 
methods  of  fitting  by  Method  1 were  attempted,  one  by  a linearisation 
technique  (lA)  C91  and  the  other  by  direct  searching  (IB). 

In  Method  3,  it  is  assumed  that  the  variance  of  y is  Independent  of 
y C13  although  there  is  no  evidence  that  this  is  so  for  the  experimental 
data  set.  For  this  reason,  the  sum  of  squares  was  weighted  as  shown  in 
Methods  4 and  5;  by  1/y  C103  or  by  l/(y+y)^.  The  latter  weighting 

factor,  which  was  used  in  another  situation,  was  claimed  nevertheless  to 
be  free  from  bias  and  applicable  to  a variety  of  enzymic  kinetic  data  Llll. 

Further  details  of  the  fitting  procedures  can  be  found  in  the 
Appendix. 


RESULTS  AND  DISCUSSION 


Description  of  Data  Sets 

Before  considering  the  results  of  the  curve  fitting  procedures,  it 
is  worthwhile  to  summarise  the  characteristics  of  each  data  set,  as 
follows.  Sets  I-III  are  simulated  data  sets  (Table  II). 

Set  I.  A = 20,  B = 19,  k = 0.05.  This  set  covers  values  of  y from 

0.26A  to  0.87A.  In  set  la,  the  sum  of  arithmetical  errors  (but 
not  the  sum  of  percentage  errors)  is  almost  zero,  while  in  set 
Ib  this  situation  is  reversed. 


Set  II.  A » 20,  B = 19,  k • 0.05.  The  values  of  the  three  constants  are 
the  same  as  in  set  I but  the  range  covered  by  the  values  of 
y (0.18A  to  0.71A)  is  not  as  great.  Sets  Ila  and  Ilb  bear  the 
same  relationship  to  each  other,  with  respect  to  the  distribution 
of  simulated  error,  as  sets  la  and  Ib  above.  Sets  I and  II  are 
not  independent  with  respect  to  distribution  of  errors. 


Set  III.  A • 20,  B » 15,  k “ 0.02.  This  set  covers  a narrower  range  still 
of  y values  (0.32A  to  0.66A).  In  set  Ilia,  neither  the  sum  of 
arithmetical  errors  nor  the  sum  of  percentage  errors  approximates 
zero,  and  both  are  negative.  In  fact,  5 of  the  8 data  points 
exhibit  a negative  error,  while  in  sets  la,  Ib,  Ila,  lib  and  IV, 
the  distribution  of  error  is  4 positive,  4 negative. 

Set  IV.  This  is  the  experimental  set  and  is  listed  in  Table  III. 

A * 27.8,  B = 21.8,  k “ 0.053.  The  values  of  y cover  the  range 
0.33A  to  0.73A.  Neither  the  sum  of  arithmetical  errors  nor  the 
sum  of  percentage  errors  approximates  zero,  and  both  are  negative. 
Set  IV  is  therefore  very  similar  to  set  Ilia  in  this  respect. 
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The  limited  number  of  data  sets  above  suffice  for  the  purpose  of 
this  paper,  which  Is  to  demonstrate  the  Inaccuracies  of  the  curve  fitting 
procedures  under  study.  Were  some  or  all  of  these  procedures  acceptable, 
many  more  data  sets  would  naturally  be  required  before  one  procedure  could 
be  recommended  in  preference  to  others. 

Results  from  Curve  Fitting 

The  curve  fitting  procedures  were  applied  to  each  of  the  6 data  sets 
under  study.  The  best-fit  values  of  the  constants  A,  B and  k for  each  of 
the  six  methods  are  given  in  Table  IV  where  they  are  compared  with  the 
"true"  values. 

Both  variations  of  Method  1 failed  to  give  reasonable  results  and 
those  which  were  obtained  were  inconsistent.  It  should  be  noted  that  the 
standard  error  S can  be  reduced  to  zero  in  this  method  by  putting 
A = B = “>  and  k = 0.  Local  minima  in  S might  be  expected,  but  this 
occurred  in  only  three  of  the  six  data  sets  using  direct  searching  (Method 
IB).  Furthermore,  the  variation  of  the  standard  error  with  the  constant  A 
in  the  vicinity  of  a local  minimum  (or  maximum)  was  erratic.  Data  set 
Ilia  actually  gave  a local  maximum  of  S for  a value  of  A (19.7)  close  to 
its  "true"  value  (20.0).  By  Method  lA,  four  of  the  six  data  sets  gave 
convergence  to  the  absolute  minimum  of  S = 0,  while  the  other  two  gave 
convergence  to  spurious  values  as  shown  in  Table  IV,  probably  because  the 
rate  of  change  in  the  corrected  values  of  the  parameters  was  slow  enough 
to  trigger  off  the  program  test  of  convergence.  The  failure  of  Method  1 
is  specially  Important  because  such  an  approach  has  been  advocated  as  part 
of  the  determination  of  the  kinetic  constants  of  a mixture  of  chemically- 
modified  and  unmodified  enzyme  C8I1.  It  should  be  noted  however  that  if  A 
is  known,  the  objection  to  this  method  no  longer  applies  as  it  then  becomes 
a simple  linear  fit  using  equation  (2).  Results  using  this  approach  are  ' ; 

shown  in  Table  VI  (see  later). 

Methods  2-5,  unlike  Method  1,  always  gave  convergence.  Superficially, 

Methods  3-5  were  satisfactory  for  data  sets  la,  Ib,  Ila  and  Ilb  in  that  the 

per  cent  deviations  of  the  best-fit  values  of  A and  B from  their  respective 

true  values  were  less  than  the  per  cent  deviation  in  some  values  of  y,  or 

nearly  so.  This  is  not  surprising  since  there  is  a uniformity  about  the 

distribution  of  error  in  data  sets  I and  II,  and  the  values  of  y span  a ] 

relatively  wide  range.  Nevertheless  the  results  are  less  satisfactory  when  j 

the  confidence  ranges  for  the  rate  constant  k are  taken  into  consideration 

(see  below  and  Table  V). 

The  results  are  certainly  unsatisfactory  in  cases  Ilia  and  IV.  In 
the  former  case.  Methods  2-5  all  gave  essentially  the  same  best-fit  value 
of  A (16.9  - 17.05)  but  this  was  lower  than  its  true  value  by  about  15%. 

Similarly  for  B and  k the  best-fit  values  differed  from  their  true  values 
by  18%  and  43%  respectively  even  though  no  "experimental"  value  of  y 
differed  from  its  true  value  by  more  than  2.33%  (a  very  acceptable  level  of 
accuracy  in  many  fields  of  science).  In  the  case  of  data  set  IV  which  is 
a real  experimental  situation,  the  maximum  "error"  in  y was  3.22%,  yet  the 
best  values  of  A,  B and  k obtained  by  any  method  were  in  error  by  6.3%, 

4.9%  and  19%  respectively. 

• s 
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Method  6 was  used  as  a trial  only  on  sets  la  and  Ilia  but  was  not 
extended  to  other  data  sets  because  of  unpromising  results  and  the  amount 
of  computing  required. 

As  a further  test  of  the  reliability  of  results  obtained  from  our 
test  data,  the  standard  deviation  of  the  fitted  value  of  k was  computed 
for  Methods  2-5  using  the  calculation  procedure  described  by  Guest  C93. 

The  values  obtained  for  each  method  were  all  much  the  same.  Table  V 
shows  the  values  obtained  by  Method  3 together  with  a 95%  confidence  range 
for  k obtained  by  assuming  that  k is  normally  distributed.  The  results 
show  how  wide  the  confidence  ranges  are  even  in  the  case  of  superficially 
accurate  values  of  k.  Where  apparently  poor  values  were  obtained  for  k, 
the  true  value  was  well  inside  the  95%  confidence  range.  It  should  be 
noted  that  the  results  in  Table  V were  obtained  on  the  assumption  of  equal 
variances  in  the  y values  so  that  the  sums  of  squares  of  residuals  could 
be  used  in  the  estimations.  The  Appendix  indicates  how  other  assumptions 
can  be  used.  Because  of  the  sensitivity  of  the  value  of  k to  experimental 
error,  the  results  reported  by  Hovanec  & Lieske  C63  (in  which  agreement  of 
the  computed  constants  A and  k with  those  observed  from  test  data  was 
claimed  to  be  within  2%)  should  not  be  taken  as  representative  of  the 
accuracy  likely  to  be  obtained  in  applications  of  this  sort.  They  used  a 
published  computer  program  L121  which  uses  our  Method  3 with  a different 
calculation  procedure  for  the  analysis  of  their  data  using  equation  (1). 

In  order  to  show  the  advantage  of  knowing  the  value  of  A,  the  best- 
fit  values  of  B and  k together  with  standard  deviations  and  95%  confidence 
ranges  for  k were  determined  for  the  six  data  sets  assuming  that  the  "true" 
values  of  A were  known  and  linear  fits  using  equation  (2)  could  be  made. 

The  results  are  shown  in  Table  VI.  They  indicate  the  better  values  obtained 
for  B and  k in  this  situation  and  the  greater  reliability  to  be  placed  on 
the  results  as  shown  by  the  tighter  confidence  ranges  for  k. 


CONCLUSION 


Our  results  suggest  that  if  an  accurate  determination  of  the  rate 
constant  in  equation  (1)  is  required,  it  is  preferable  to  concentrate  on 
determining  the  constant  A experimentally,  if  necessary  by  an  indirect 
method.  An  alternative,  which  would  probably  be  more  time-consuming,  is 
to  make  a thorough  investigation  of  the  nature  of  the  experimental  error 
in  the  system  under  study  C13I1;  however  knowing  this,  the  results  of 
curve  fitting  procedures  may  still  not  be  sufficiently  accurate.  This 
conclusion  is  important  when  it  is  realised  that  equation  (1)  is  not  con- 
fined to  enzyme  kinetics,  but  describes  many  processes  in  other  fields  of 
science  C143  and  engineering  C123.  Curve  fitting  procedures  have  been 
adopted  in  these  references  C 12, 14 3 and  also,  in  the  field  of  enzyme 
reactivation,  by  Lieske  et  al.  C153. 
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Ilia  y=y  + (E3)  6.286  7.817  8.786  10.177  10.884  11.902  12.392  13.090 
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used. 


The  following  are  brief  descriptions  of  the  curve  fitting  methods 


Method  1 

This  method  alms  to  minimise  the  standard  error  of  the  least-squares 
line  through  the  n data  points  lIln(A-yp  ,t  .1.  The  minimum  of  S is  sought 
by  determining  appropriate  values  of  A,  B and  k where 


as 

= Cln(A-y^)-lnB  + kt^3^. 


i«l 

As  stated  In  the  text,  S can  be  reduced  to  zero  by  putting  k=0,  B=A  and 
then  allowing  A to  Increase  Indefinitely.  As  this  Is  an  unacceptable 
result  for  our  purposes,  a local  minimum  occurring  elsewhere  Is  sought. 
Two  methods  were  used  In  an  attempt  to  achieve  this.  The  first  method 
(Method  lA)  is  a linearisation  technique  CID.  This  starts  with  approxi- 
mate values  A , B and  k of  the  respective  parameters  and  then  tries  to 
find  small  (assumed)  corrections  6A,  5B  and  6k.  In  this  experiment,  the 
Initial  values  used  were  the  "true"  values  of  the  parameters.  Following 
Guest's  procedure  C13,  the  sum  of  squares  to  be  minimised  becomes 


where  6A,  6B,  6k  are  the  quantities  to  be  determined  and 


yi^  = ln(A^-y^)-lnB^4k^t^  . . 

i 

j 

Then  Improved  values  of  A,  B,  k are  determined  by  adding  the  corrections 
to  the  starting  values.  If  more  accuracy  Is  required,  the  procedure  can 
be  Iterated  as  often  as  desired  using  the  Improved  values  as  new  Initial  l 

values . j 

In  the  other  method  (Method  IB)  the  minimum  of  S was  sought  by  a | 

direct  search.  In  which  values  of  A were  the  basis  of  search,  corresponding  ' 

values  of  B and  k being  obtained  by  the  usual  linear  fitting  technique.  \ 

Because  of  the  failure  to  achieve  useful  results.  Methods  lA  and  IB  were 
abandoned  as  being  of  little  use  for  this  type  of  problem.  j 


= ■ s ) 

1-1  ' 

This  problem  Is  discussed  fully  by  Guest  HID  who  uses  the  linearisation 
process  outlined  under  Method  lA.  Overall  the  results  using  this  method 
were  the  best  obtained.  It  gave  the  best  results  for  data  sets  Ia,b  and 
Ila.  It  was  not  far  from  best  for  data  sets  Ilb  and  Ilia. 

Method  4 

Similar  to  Method  3 but  using  weights  where 

y^  - A-Be 

Not  as  good  overall  as  Method  3 for  the  data  sets  used. 
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Because  of  the  awlcward  nature  of  S,  direct  search  methods  seem  to  be 
necessary.  The  amount  of  computing  Is  large  but  can  be  reduced  by  the 
following  procedure. 

Consider  B and  k as  fixed  quantities  and  calculate 

-kti 

= y^+Be  , 1 = 1,2,  , n . 


Now  rearrange  the  values  of  as  where 


< C2  < 


< C 

n 


Then  for  the  particular  values  of  B and  k,  S Is  minimum  for  A equal  to  the 
median  of  the  C . The  supreme  minimum  Is  then  found  by  searching  on  B and 
k.  The  method  was  used  only  on  data  sets  la  and  111a,  but  as  results  were 
not  worth  the  amount  of  computing  Involved,  the  method  was  not  tried  with 
the  other  data  sets. 

2.  ERROR  ESTIMATES 

The  usual  method  of  obtaining  an  estimate  of  variance  In  least-squares 
problems  Involves  computing  the  residual  sum  of  squares  using  the  fitted 
parameter  variables.  This  procedure  Is  described  by  Guest  and  Is  applicable 
to  all  of  Methods  1-5  described  above.  The  variances  of  the  parameter 
estimates  can  then  be  determined  and  hence  confidence  ranges  can  be 
calculated  If  assumptions  about  the  statistical  distributions  of  the 
estimates  are  made,  e.g.  that  they  are  normally  distributed.  This  procedure 
however  Involves  the  assumption  that  the  variance  (and  If  applicable  the 
forms  of  distribution)  of  the  y values  Is  the  same  throughout.  It  Is  quite 
possible  for  this  not  to  be  so.  In  this  case  the  following  procedure  may 
be  useful.  It  refers  to  Method  3 but  can  easily  be  adapted  for  other 
methods. 
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The  normal  equations  have  the  solution 


where 


6A  - (QX-PY+NZ)/D, 

6B  ■ (PX-MY+LZ)/D, 

6k  - (NX-LY+KZ)/B  D 
o 

D - nQ-SP+TN, 

K - nU-S^,  L-nV-ST,  M - nW-T^, 
N - SV-TU,  P - SW-TV,  Q - UW-V^ 


and 


_fL  -k  t. 

n 

-V  ^ _5L 

2^'  ■ ^ 

-L 

^1® 

- z 

1-1 

1-1 

1-1 

-2k  t 

n 

, 2 
h ® 

. 0*^1 

w » 

z 

1-1 

1-1 

-2k  t , 
o 1 


1-1 


2«yi,  Y - 2^e  ° Vi.  Z - ° Vi  , 


-k  t. 


1-1 


1-1 


where  A , B , k are  the  Initial  estimates  of  A,  B,  k respectively  and 
6A,  6B,  6k  are  °the  corrections  to  be  added.  The  6y^  are  defined  by 


6y  - y.-A  +B  e 
•'I  •'loo 


-k  t, 
o 1 


The  solution  can  be  expressed  In  the  form 


n 


6A 


^A^6y^,  6B  - 5^B^6y^.  6k  - k^6y^ 


1-1 


1-1 


1-1 


(1) 
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where  the  coefficients  In  the  summations  are  given  by 


DAi 

-k  t 

- Q-(P-Nt^)e  ° ^ 

DBi 

-k  t . 

- P-(M-Lt^)e  ° ^ 

-k  t. 

- N-(L-Ktj^)e  ° ^ 

Any  Information  or  assumptions  about  the  variability  of  the  y^,  which  is  the 
same  as  that  of  the  6y. , can  be  used  In  the  usual  ways  to  obtain  correspond- 
ing variabilities  of  the  parameter  estimates.  For  example.  If  the  variances 
of  the  y^  are  the  variance  of  the  estimate  of  k will  be  estimated 

In  the  usual  way  by 


n 


E 

1-1 


2 

1 


k 


2 

1 


and  similarly  for  the  other  parameters.  Systematic  errors  can  also  be  dealt 
with  If  required. 

As  a numerical  example,  data  set  Ilia  gives  for  the  5k  expansion, 
using  the  "true"  values  of  A,  B,  k for  A , B , k 


100  5k  - - 1.721  6y^  + 0.015  5y2  + 0.983  &y^  + 1.312  5y^  (2) 

+ 1.116  5y^  + 0.492  5yg  - 0.477  5y^  - 1.720  5yg 

and  similar  expressions  for  5A,  5B. 

From  (2),  results  such  as  the  following  can  be  derived. 

(a)  Substituting  the  actual  values  of  the  5y.  for  this  data  set  gives  the 

value  0.00862  for  5k,  the  Increment  to  be  added  to  k on  the  first  Iteration 
of  method  3,  giving  k » 0.0286.  ° 

(b)  If  the  greatest  possible  error  In  the  y^  is  ± 0.25  and  the  errors  are 
cumulative,  then 


5k 


< 0.25 


1-1 


giving  a possible  range  for  k of  0.0286  ± 0.0196. 
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(c)  If  It  Is  assumed  that  the  standard  deviations  of  the  are  the  same 

and  that  the  maximum  error  Is  ± 0.25  and  that  this  represents  the  95% 


confidence  limits  for  the  y.  values,  then  the  variance  of  6k  and  hence 


of  the  k estimate  Is  given  oy 


n 


0.25 


1.96 


0.0000168 


1=1 


so  that  the  95%  confidence  range  for  k Is  0.0286  ± 0.0105. 

(d)  If  It  Is  assumed  that  the  standard  deviations  of  the  y^  are 
proportional  to  the  true  values  of  the  y^^  and  that  the  95%  confidence  range 
for  each  y^  Is  ± 2.5%  of  the  true  value,  then 


0,025 


1.96 


k,^y.^  = 0.0000180 


^*"1  ^1 


so  that  the  95%  confidence  range  for  k Is  0.0286  ± 0.0109.  In  an 


experimental  situation,  the  "true  values  of  the  y^  would  be  unknown,  but 
In  this  case  the  experimental  values  could  be  used  without  changing  the 
results  significantly. 
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